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Abstract 

We construct the fundamental solution or Green function for a divergence 
form elliptic system in two dimensions with bounded and measurable coef- 
ficients. We consider the elliptic system in a Lipschitz domain with mixed 
boundary conditions. Thus we specify Dirichlet data on part of the boundary 
and Neumann data on the remainder of the boundary. We require a corkscrew 
or non-tangential accessibility condition on the set where we specify Dirichlet 
boundary conditions. Our proof proceeds by defining a variant of the space 
BMO{Q) that is adapted to the boundary conditions and showing that the 
solution exists in this space. We also give a construction of the Green func- 
tion with Neumann boundary conditions and the fundamental solution in the 
plane. 
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1 Introduction 



We consider a weak formulation of the mixed problem for a second-order elliptic 
operator in a bounded, connected, and open set f2 in R 2 . To state the mixed 
problem, we fix a decomposition of the boundary d£l = D U N with D n N = 
0. We let L be an elliptic operator in divergence form with bounded and 
measurable coefficients and we consider the boundary value problem 

Lu = f in Q 

u = on D (1.1) 

$ = f N , on N. 

where d/du is the natural co-normal derivative associated with the operator L. 
We will require that be a Lipschitz domain and that the set D satisfy a cork- 
screw condition (or non-tangential accessibility condition) as in [25] . We will 
give a precise formulation of the mixed problem in section [2j We emphasize that 
our results apply to elliptic systems where L acts on vector-valued functions 
as well as equations. Our goal is to give a construction of the matrix-valued 
Green function for this boundary condition and show that the Green function 
satisfies the estimates 

\G(x,y)\ < 
\G(x,y)-G(x,z)\ < 

\G(x,y)\ < 

where d is the diameter of and the constant C and 7 > depend only 
on assumptions on the operator, the domain, and the decomposition of the 
boundary. 

Given the logarithmic singularity of the Green function, it is natural to look 
for a Green function in the space of functions of bounded mean oscillation. We 
fix a set D and define the space BMOd{Q) which consists of functions in 
BMO{Q) and which vanish in an appropriate sense as we approach D. We 
will construct the Green function in this space and then obtain the pointwise 
estimates. 

There is a great deal of literature on the existence of Green functions and we 
do not try to summarize it all here. Littman, Stampacchia, and Weinberger [18] 
establish that the Green function has a logarithmic singularity in two dimen- 
sions. Griiter and Widman p3] give a nice construction of the Green function 
for the Dirichlet problem in dimensions three and higher and treat operators 
that are not self-adjoint. Kenig and Pipher [T7] give a construction of the Neu- 
mann function or the Green function for the Neumann problem when n > 3. 
The existence of a global fundamental solution for an elliptic operator in the 
plane was established by Kenig and Ni |16j . Using their results, Chanillo and 
Li [3] show that the global fundamental solution lies in BMO(R 2 ). 
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The results of the previous paragraph are for single equations. More re- 
cently, there has been interest in constructing Green functions for elliptic sys- 
tems. In this case, we are only able to obtain upper bounds. The work of 
Auscher and collaborators [T] establishes the existence of a Green function for 
operators with complex coefficients (which may be viewed as a system of equa- 
tions with real coefficients). This work was motivated in part by their interest 
in the Kato problem. Dolzmann and Midler [7j construct the Green function 
in a C 1 -domain for an elliptic operator with continuous coefficients. Dong and 
Kim [8j establish the existence of a Green function for elliptic systems which 
are not assumed to be self-adjoint under the hypothesis that solutions of the 
operator are bounded or Holder continuous. Their argument gives Gaussian up- 
per bounds for the parabolic Green function and then integrates the parabolic 
Green function in time giving the Green function for the elliptic problem. Re- 
cently, Choi and Kim [5J have given a construction of the Neumann function in 
dimensions three and higher under the assumption that solutions of the operator 
satisfy local Holder continuity and boundedness estimates. In two dimensions, 
the necessary results on the Holder continuity of solutions dates back to Mor- 
rey [21] . A paper of Calanchi, Rodino, and Tri [2] observes that the Green 
function for the Dirichlet problem for the Laplacian lies in BMO for planar 
domains. Their proof relies on the maximum principle which is not available 
to treat the systems considered here. In addition to our main goal of treating 
the Green function for the mixed problem, the present paper complements the 
work of Choi and Kim by providing a construction of the Green function for the 
Neumann problem in two dimensions. Recent work of D. Mitrea and I. Mitrea 
[20| gives the existence of the Dirichlet Green function for various constant 
coefficient operators in Lipschitz domains. 

The present paper was motivated by an effort to construct the Green func- 
tion for elliptic systems under mixed boundary conditions. Taylor, Ott, and 
Brown |25| give a construction of the Green function for a class of mixed prob- 
lems for the Laplacian in a Lipschitz domain in dimensions two and higher. 
Their argument begins with the free space fundamental solution and uses a 
reflection to first construct a Green function for Neumann boundary conditions 
and finally corrects the boundary conditions to give a Green function for mixed 
boundary conditions. Below, we provide a new proof of their results in two di- 
mensions and establish the existence of a Green function with mixed boundary 
conditions for a large class of elliptic systems in two dimensions. We expect 
that the Green function constructed below will be useful in extending Taylor, 
Ott, and Brown's results on the L p -mixed problem for the Laplacian to systems 
in two dimensions. It is an interesting, open question to study the mixed prob- 
lem for elliptic systems in dimensions three and higher. The main difficulty in 
carrying out this extension is the lack of estimates for the Green function. 

The construction of the Green function for an operator L is closely related 
to local, scale-invariant estimates for solutions of L. In the case of a single equa- 
tion, the Holder continuity of solutions to the mixed problem was established 
by Stampacchia [23] using the method of De Giorgi. Stampacchia considers 
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the mixed problem in domains which locally are equivalent by a bi-Lipschitz 
transformation to a mixed problem in a half-space with the boundary between 
D and ./V a hyperplane. A more recent work of Groger [13] uses the method of 
Meyers [19] to show that solutions of the mixed problem with nice data satisfy 
Vii € L 2+e (Q,). Groger 's assumptions on the domain and the decomposition 
of the boundary are similar to those of Stampacchia. When considering the 
more general decompositions of the boundary introduced by Taylor, Ott, and 
Brown and considered in this paper, it seems to be simpler to use the method 
of reverse Holder inequalities as in Gehring [9] and Giaquinta and Modica |11] 
to obtain that Vu G L 2+e (fi). 

Our formulation of the mixed problem allows for the extreme cases where 
D = or D = dQ which give the Neumann problem and the Dirichlet problem, 
respectively. The Green function for the Dirichlet problem is treated alongside 
the mixed problem in section [U The properties of the Green function are given 
in Theorem 14.11 The changes needed for the Neumann problem appear in 
section [5] and the properties of the Green function for the Neumann problem 
are given in Theorem 15.141 In Section [6l we provide a proof of the existence of 
a fundamental solution in the plane. The construction of a global fundamental 
solution for an equation was perhaps known, however the detailed construction 
appears to have first been written down in a paper of Kenig and Ni from 1985 
[16j . We provide a different construction of the fundamental solution which 
also applies to systems. The existence of a fundamental solution in the plane 
and the properties of this fundamental solution are given in Theorem 16. 71 The 
Lame system (with variable Lame parameters) provides a family of examples 
to illustrate the use of our results. We describe this system in section [2] and 
show that it satisfies the hypotheses of our main results. 

2 Preliminary material and the weak for- 
mulation of the mixed problem 

We will consider the Green function for boundary value problems in a bounded 
Lipschitz domain and we begin by giving the definition of these domains. 
Given a constant M > 0, x 6 <9Q, and r > 0, we let Z r {x) = {y : \y± — x±\ < 
r, | j/2 — ^2 1 < (3M + l)r}. We say that Z r (x) is a coordinate rectangle for $7 if 
there is a Lipschitz function : R — > R so that 



We assume that the coordinate system used in this coordinate rectangle is 
a rotation of the standard coordinate system. We say that f2 is a Lipschitz 
domain if f2 is a bounded connected open set and for each x € d£l, we may find 
a coordinate rectangle centered at x. Since the boundary is compact, we may 
cover <9f2 by a finite collection of coordinate rectangles {Z r . (xi)}^ =l so that each 



o n z r (x) 
dnnz r (x) 



{y-V2> 4>(yi)} n Z r (x) 
{y-V2 = 4>(yi)} n z r (x) 
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ZiQOri{xi) is also a coordinate rectangle. We set ro = minjrj : i = 1, . . . ,N} 
and we will use this value as a characteristic length of the domain when we 
state scale- invariant estimates. 

For a; E ft and p € (0, 4ro), we define local domains ftp(x) which will play the 
role of disks in our work. In the case that dist(x, <9ft) > p, we let £l p (x) = B p (x), 
the disk centered at x with radius p. In the case that dist(x, <9ft) < p, we pick a 
coordinate rectangle containing x = (xi,X2) and using the coordinates for this 
rectangle, we let x = (xi,(f>(xi)) be obtained by projecting onto the boundary 
and define £l p (x) = Z p (x) n ft. 

For x E c?$7 and p £ (0, ro), we define a boundary interval A p (x) by 
A p (x) = Z p [x) n <9ft. The domains Q p (x) are star-shaped Lipschitz domains 
and the boundary intervals A. p (x) are connected. These properties are helpful 
in establishing the Poincare and Korn inequalities introduced below. For this 
reason, we prefer them to B p (x) n ft and B p (x) n <9ft used in [231 EE], for ex- 
ample. There is a price to pay as given x and p we will have several choices 
for A p (x) and ftp (a;). Our results will hold for any such choice provided that 
when several of these objects occur in an estimate we use the same coordinate 
rectangle to define each of them. 

We let D C <9ft be the set where we specify the Dirichlet data and then 
we put N = <9ft \ D. We will require that D satisfy the interior corkscrew 
condition. This means that for each x E 3D (where the boundary is taken with 
respect to 9ft) and r £ (0, ro), we may find x r € D so that 

\x-x r \ < r and dist(x r , JV) > M _1 r. (2.1) 

As our boundary is locally a Lipschitz graph, this easily implies that for each 
x G D and r £ (0,ro), we have that a(A r (x) C\ D) > cr (see [25], Section 1] for 
details) . Note that if D is not empty, then the corkscrew condition implies that 
the interior of D is not empty. 

Now we turn to the precise description of the boundary value problem. 
As we are considering an operator acting on vector- valued functions, most of 
the functions we consider will take values in R m for some m. We do not 
explicitly denote the target in our notation for function spaces. However, we 
emphasize that in the definitions below, all of the functions will be vector- 
valued. We let W 1,p (£l) denote the homogeneous Sobolev space of functions 
having one derivative in L p (£l) with the norm = II ^ u \\lp(q)- Given a 

set D C 30, we let W^ p (Cl) denote the subspace obtained by taking the closure 
in W 1,p {£i) of the smooth functions in Cl which vanish in a neighborhood of D. 
Thus, the elements in W^ P (Q) are functions which in some sense vanish on D. 
When D is not empty and satisfies the corkscrew condition and 1 < p < oo, we 
have that if u £ W^ p (n), then 

r Q 1/p ( [ \u\ p da) + r 2/p ( [ \u\ p dy) ^ < C ( [ \Vu\ 2 dy) ^ . (2.2) 
\Jan J \Jn J \Jn J 

for a constant C which depends on M, p, and 0. We will let W' 1 ' 1 (ft) be 
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the dual of Wjf(tl) when 1 < t < oo. The space W]j ' (dO.) is defined 
to be the image of W^ 2 (fl) under the trace map and then W D ^ 2 ' 2 (d££) the 
dual of Wp /2 ' 2 (dQ). We let (-,-} : W£ 1,2 (fi) x Wj ,2 (n) -»• R and ( v )an : 
Wn 1/2 ' 2 {dn) x W# 2,2 (0fi) ^ R be the pairings of duality. 

The operator L will act on vector-valued functions in the plane. Formally, 
L is given by 

N°=EE^fe a = l,...,m. (2.3) 

i,j=i p=i 1 i 

The coefficients a£p are assumed to be real, bounded, and measurable functions 

max {||«a/3llL°°(o) : h3 = 1,2, a, (3 = 1, ... , to} < M. (2.4) 

We do not assume a symmetry condition on the coefficients and thus the op- 
erator L will not be self-adjoint. We will define Lu as an element of W~^{Vt) 
via the quadratic form 



dvP dcj) a 



Here and throughout this paper, we follow the convention that we sum on 
repeated indices. Since the coefficients are bounded, for some constant C = 
C(M, to) we have 

\A(u^)\ <C||«|| v j rl , 2(n) ||^|| T ^ 1 , a(n) . (2.5) 



We assume the following local ellipticity condition. If x € fi, p E (0, ro), and 
u € Wft' 2 (0,), then for all constant vectors c € R m , we have 

M- 1 [ \Vu\ 2 dy<[ a%^^ + p- 2 \u-c\ 2 dy. (2.6) 
Jnjx) Jn B (x) P9 Vj % 



du@ du a 

In addition, we assume a global coercivity condition on D and the form A 



M- 1 [ \Vu\ 2 dy <A(u,u), ueW^ 2 (n). (2.7) 



The conditions (|2.6p and (|2.7p are immediate if D has non-empty interior 
and the coefficients satisfy the strong ellipticity condition 

<>';!a;c; > ciei 2 , e e R 2m 

At the end of this section, we will show that they also hold for the Lame system 
when the set D is non-empty and satisfies the corkscrew condition. 
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We are now ready to give the precise formulation of (jl.ip . Given / in 

12 "1 /2 2 

WV' (°) and fN G * (90) we say that it is a weafc solution of the mixed 

problem if we have 

0) = -(/, 0) + 0)an, e V^' 2 (ft) 

With our continuity ()2.5j) and global coercivity (|2.7[) assumptions, the existence 
and uniqueness of solutions is an immediate consequence of the Lax-Milgram 
theorem. 

We also consider the adjoint problem for the operator L* whose coefficients 

ji 

mixed problem for L* 

L*w = g, in Q 

^=g N , onN 
w = 0, on D 



are obtained by replacing a^o by a J o Q . We say that w is a weak solution of the 



if we have 



A((j),w) = -(g,0) + (g N ,0) m , ^eW 1 /^) 



w 6 W£ r(o). 

Given a locally integrable function on CI, x E Cl, a subset D C 917, and 
p G (0, ro), we define 

if dist(n p (x),D) = 
u(y)dy, if dist(O p (a?), Z>) > 

Here, we are using f E f dy = J B / <iy to denote the average of E. We will 

define the space BMOo{Q) to be the collection of integrable functions wonU 
for which the norm 

IMI*,D = sup{-f | it - u X)P \ dy : x G Cl, p € (0, r )} 

is finite. The supremum is taken over all local domains. We will be primarily 
interested in the case when D ^ and then it is easy to see that the above 
expression is a norm. In the case where O is R 2 or D = 0, then the elements 
of the space will be equivalence classes of functions which differ by a constant. 

We note that if D 1 C D 2 , then we have that BMO Dl {9) D BMO D , 2 (Sl) 
and in particular we have BMOd{Q) C BMO${VL). Thus, we obtain the John- 
Nirenberg inequality 

\{y £ ttp(x) : \u(y) - u x>p \ > A}[ < Cp 2 exp(-A/||it||* )D )}. 

The standard proof for BMO(R n ) as found in the lecture notes of Journe 
[15j Chapter 3], for example, extends easily to the space BMOd(&)- As a 
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consequence, it follows that if D ^ 0, then BMO D {9) C LP(Q), for 1 < p < 
oo. In the case where f2 = R 2 , we can show that any representative of a 
BMO(R 2 ) function lies in n p<00 L^ 0C (R 2 ) and in the case of a bounded domain 
any representative of a function in BMO${Q) also lies in n p<00 L p (£l). 

We define atoms for Q and D and then the Hardy space iJp(O). We say 
that a bounded measurable function a is an atom for £1 and D if a is supported 
in one of the local domains Cl p (x) and satisfies 

Nlr<»(n) < 1 /\ n p( x )\ 
a x , P = 0. 

A function / is in the atomic Hardy space Hp(Q) if there is a sequence of atoms 
{cii}?^ and a sequence of real numbers {Aj} € i 1 so that / = Yl*S=i ^i a i- We 
define a norm on this space by 



JJ* (O) 



inf ^ | Aj 



where the infinum is taken over all representations of /. 
It will be useful to observe that the expression 

sup{ / a a u a dy : a is an atom for Q and D} (2-10) 
Jn 

gives an equivalent norm on BA10e>(£1). This proposition may be found in 
Journe [El Chapter 3] for BMOiW 1 ). The extension to BM0 D {9) is straight- 
forward. We will use the characterization of the BMOo{Q) norm in (|2.10|) to 
show that the Green function lies in BMOd(Q)- In fact, BMOo{Q) may be 
identified with the dual of the atomic Hardy space Hpifl). See Journe [To] . 
Coifman and Weiss [6] for the case when D = 0. Chang [3] treats the extreme 
cases where D = or D = <9$7. The extension to general D is not difficult. 

We recall several Poincare and Sobolev inequalities that will be needed in 
the argument below. The first inequality is a scale invariant Sobolev-Poincare 
inequality. If u G W/ifL), 1 < p < 2 and 1/q = 1/p - 1/2, and fi p (x) is one 
of the local domains defined above, then we may find a constant C depending 
only on M and p so that 

-u x , P \ q dy) <c( f \Vu\Pdy\ . (2.11) 

If u X)P = fa ^udy, the inequality holds with f2 p (x) as the domain of inte- 
gration on the right-hand side. However, if dist(f2 p (x), D) = 0, then we need 
to expand Q p (x) and use the corkscrew condition (12. ip in order to conclude 
that u vanishes on a large enough set to obtain the inequality (|2.1ip . See [231 
Section 3] for more details. A useful consequence of the inequality (12. llf) is 
the Poincare inequality for 1 < p < oo, 

\ Vp / . \ Vp 

u x JPdy) <Cp( \Vufdy) (2.12) 

/ \Jn 2p (x) J 





Note that the Poincare inequality ()2. 12j) with p = 2 immediately implies the 
embedding Wrf (0) C BMOd(&)- The corkscrew condition is needed to estab- 
lish this embedding. If 2 < p < oo, we have the following version of Morrey's 
inequality 



u — u 



*,plk«(n,(*)) < CP 1_2/P / |Vnr . (2.13) 




This inequality gives the embedding of the Sobolev space Wjf(n) into space 
of Holder continuous functions with exponent 1 — 2/t when t > 2. Finally, we 
give an estimate at the boundary. Let £t p (x) be one of our local domains and 
suppose that l<p<oo, l<g<2 and 1/q = l/(2p) + 1/2. There exists a 
constant C = C(p,M) so that for u G W^ q (Q,2p(x)) , we have 

(f \u-u x JPdy] <c(f \Vu\"dy) . (2.14) 
\Ja p (x) J \Jn 2p ( x ) J 

To end this section, we recall the Lame operator and show that the form 
for this operator with the mixed boundary condition satisfies the ellipticity 
condition (12. 6h and coercivity condition (]2.7h . Given two real- valued functions 
[i and A, the Lame operator is the operator with coefficients 

a ap( X ) = M X )( ( %<W + Sj a 5ij3) + X(x)5 ia 5 jl 3 
We will use e(u) to denote the strain or the symmetric part of the gradient, 

, . , If du l du° 



2 9xj 

and then cr(u) will denote the stress tensor which is given by 
or more compactly by 

o~{u) = 2 [i e(u) + A I div u 

where J is the 2x2 identity matrix. The functions \i and A are called the Lame 
parameters and are related to the elastic properties of the material. We require 
that A and \x are bounded, measurable functions and that \i{x) — A _ (x) > c > 
where A - denotes the negative part of the function A. With this assumption, 
we have the pointwise lower bound 

a^>c\e{u)\\ (2.15) 
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We recall the following version of Korn's second inequality. Let £l p (x) be 
one of our local domains with star-center x*. If u is in W ' 2 (Q p (x)) and c is 
any constant vector in R 2 , then 

/ \Vu\ 2 dy<C([ \e{u)\ 2 dy+\( \u - c\ 2 dy) (2.16) 

Jn p (x) Jn p (x) P Jb p/2 (x*) 

The constant depends only on M. This may be established using the argument 
given in the monograph of Oleinik, Shamaev, and Yosifian \22\ Theorem 2.10]. 
It is clear that the ellipticity condition (12. 6p follows from the pointwise bound 
(|2.15p and (|2.16p . If D is non-empty (and then the corkscrew condition implies 
D has non-empty interior), the coercivity condition (|2.7|) may also be found in 
Oleinik, Shamaev, and Yosifian [22\ Theorem 2.7], however the standard proof 
of this inequality seems to be by contradiction and thus we cannot say anything 
about how the constant depends on the domain Q and the boundary set D. 

All of our quantitative assumptions have been in terms of the constant 
M. The results below will be of two types. Many of the results will be local 
estimates which hold on scales p S (0, tq). In these local estimates, the constant 
will depend only on m, M, and any L p -indices that appear in the estimate. The 
remaining estimates will depend on global properties of the domains such as 
the collection of coordinate cylinders which cover the boundary or the constant 
in the Poincare inequality (|2.2p . However, the constant may be chosen to 
be uniform under small changes in the Lipschitz functions which define the 
boundary. One exception is the constant in the Korn inequality which gives the 
coercivity condition (j2.7|) for the Lame system. As noted above, the standard 
proof is by contradiction and gives no information about the dependence of 
the constant on the domain. In the study of the Neumann problem the global 
estimates for solutions will depend on the estimate in our existence theorem, 
Theorem 15.41 The proof of this result uses the Fredholm theory and thus we 
have no information about the behavior of the operators used to solve the 
Neumann problem. 

3 Estimates for solutions of the mixed prob- 
lem 

In this section, we prove estimates for solutions of the mixed problem. We 
will use the reverse Holder technique of Giaquinta and Modica [11] . Before 
beginning the main argument, we introduce two auxiliary functions which will 
arise when we prove Caccioppoli-type inequalities for solutions of the mixed 
problem (|2.8p when the data / or is not zero. 

In the following definitions, we assume that < r < ro so that we can define 
fip(x) and A p (x) for < p < r. For < a < 2, we let I a ,rf(x) denote the 
maximal fractional integral given by 

I a ,rf(x)= sup p a ~ 2 / \f{y)\dy. 
0<p<r Jn p (x) 
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Note that we have several choices for Q p (x) when x lies in several coordinate 
rectangles. In defining I ar , we take the maximum over the choices for £l p (x) 
arising from our finite cover of the boundary by coordinate rectangles. We have 
that p a ~ 2 Xn p (x)(y) < C(M)\x-y\ a ' 2 and it follows that I a , r f(x) < CR a \f\{x) 
where R a is the standard Riesz potential. Thus, the Hardy-Littlewood-Sobolev 
theorem gives us that with q and p related by 1/q = l/p + a/2 and 2/(2 — a) < 
p < oo, we have 

UWIUp(n) <C(p,M)||/|| w(n) . (3.1) 
Next, if is a function on <9f2, we define 

PrfN(x) = sup{4 |_fjv| da : dist(x, dQ) < p < r} 

with the convention that the supremum of the empty set is zero. Recall that 
if x is in a coordinate rectangle Z r (y), and x = (xi,<p(xi) + te2), we defined 
x = (xi, 4>{x\)). We also note that if dist(x, dVL) < p, then A p (x) = d£lnd£l p (x) 
and if dist(x, <9S1) > p, then dQ,ndd p (x) = 0. Our Neumann data, /n is initially 
only defined on N. In the definition of P r , we assume that fjy has been extended 
to d£l by setting /at = outside N. Following the argument in [231 Section 3], 
for 1 < p < oo, we may find a constant C = C(M,p,Q) so that we have the 
estimate 

r \ l/2p / p \ 1/p 

/ |P,/| 2p dy <C / l/^drr ) (3.2) 

and the corresponding estimate for p = oo is trivial. 

We are now ready to begin our estimates. We begin with a simple energy 
estimate. 

Proposition 3.3 If f lies in the Hardy space Hj-,(Q,) and f^ = 0, then the 
weak mixed problem \2. 8\) has a unique solution u which satisfies 



[ \Vu\ 2 dy<C\\f\\ 2 Hl 
Jn D 



(ay 



Proof. Let a be an atom for the Hardy space Hp(Q). Using the Poincare 
inequality (I2.12j) . it is easy to see that the map u — > Jq a a u a dy lies in the dual 
of Wp 2 (Q.). As a consequence, the solution of the mixed problem (|2.8p with 
/ = a and /at = will satisfy \ Vu\ 2 dy < C . The result for a general element 
of Hjj(Q,) follows easily from the result for an atom. I 

The following theorem gives our main estimate for solutions of the mixed 
boundary value problem. We note that the estimates of this theorem hold for 
the extreme cases where D = or d£l. In the arguments below, we will use the 
standard notation 1/p' = 1 — 1/p to denote the conjugate exponent. 

Theorem 3.4 Let Q be a Lipschitz domain, let D C dVt be a (possibly empty) 
set that satisfies the corkscrew condition \2. and suppose that L is an elliptic 
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operator which satisfies the conditions \2.J$ and \2. 6\) . Let u be a weak solution 
of the mixed problem \2. 8\) and assume that f E L P (S1) and fjy E L PN (N) are 
functions. We may find a to > 2 so that with t G [2, to), p G (1,2) and 
Pn G (l,oo), and p G (0,ro) we have the estimate 



|Vit|* dy 



< C 



1/2 



\Vu\ 2 dy 



+ 



(/ 



1/f 



P2 P (\f N \ PN ) t/pN dy 



The exponent to depends on M, p, andp^ and the constant C = C(M,t,p,pN). 

Proof. We let u be a weak solution of the mixed problem (I2.8P with / G L P (Q) 
and /at G L Pn (JV) for some p > 1 and p^r > 1- Define g by 1/g = 3/2 — 1/p 
and gA? by 1/gA? = 1 — 1/(2pn) and let s = max(q, g^r) < 2. We fix x G fi, 
< p < ro and we claim the Caccioppoli inequality 



/ 

Jn, 



\Vu\ 2 dy 



< c 



\Jn 4p (x) 



2/s 



|V< dy + I p , 2p (\f\ P ) 2/p (x) + P2 P (\fN\ PN ) 2/m (x) 



(3.5) 



The estimate of the theorem follows immediately from (|3.5p and the reverse 
Holder argument of Giaquinta and Modica (see [10} Theorem 1.2, p. 122]). To 
establish ()3.5|) . we fix x G Q and p G (0, ro). We let rj be a cutoff function 
which is one on Q p (x), zero on Q \ Q,2 P (x), and satisfies |Vr/| < C/p. We let 
£ = J* n2 ^ |u — u X: 2 P \i]\Vr]\ |Vu| + p~ 2 \u — u Xt 2 P \ 2 dy. As our first step, we will 



establish the inequality 



\ [ ri 2 \Vu\ 2 dy<C \\e + 
P Jn 2p (x) \ P 2 



2/s 



|Vnj s dy 



+ 



,2-p 



f2 2p (x) 



\ 2/p 

l/l p dy 



Q 4p (a;) 

7" 



9f2nf22p(^) 



i/ivr 




(3.6) 



If we use (|2.1ip and the Cauchy inequality with e's, the estimate (|3.6p implies 
(13. 5p . Thus we turn to the proof of (13. 6p . 

To establish (|3.6p . observe that thanks to our definition of u x „, we have 



v = rj 2 (u — u Xi 2 P ) G W^,' 2 (r2) whenever u G W^(fi). Thus, from the product 
rule, our ellipticity assumption and the weak formulation of the mixed problem, 



1.2, 
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(|2.8p . we obtain 

/ rf\Vu\ 2 dy< [ \V(ri(u-u x , 2p ))\ 2 dy + CE 
Jn Jn 

< C(A(r](u - u Xt 2p),rj(u - u x ^ P )) + E) 

< C(A(u,v) + E) 

< C( / da - [ f a v a dy + E) 

Jan Jn 



(3.7) 



We claim that if 1 < p < 2, we have 



-4 / f a v a dy 

P Jn 



< 



(/ 

yn 2p {x) 



\ 2/9 

|Vu| 9 dy] +C 



(-5=- /" l/r^yl • (3- 

\P 2 p Jn 2p (x) J 



We also will need that if 1 < pn < oo, then 
Jn 



< 



/ iVti^dj/l +C(/ |/ivP<H • (3.9) 

Jn 2f ,(x) I \J dnndn 2p (x) 



It is easy to see that (|3.6[) follows from (|3.7H3.9|) and thus it remains to prove 

(prgp and dnsj). 

To establish (|3.8|) . we use Holder's inequality with 1 < p < 2 and then the 
Sobolev-Poincare inequality (|2.1ip to obtain that with 1/q = 3/2 — 1/p 



\ [ f a v a dy 
P 2 Jn 



1 

< T2 



i/p 



< 



P \Jn 2p (x) 
C 



\f\ p dy 



i/p> 



\u - u x ,2 P \ p dy 



p 2/ q -l+2/ P I J 



n 2p (x) 
i/p 



\f\ p dy 



1/9 



|Vu| 9 dy 



n 4p (x) 



Now (I3.8p follows easily. 

To establish (|3.9p . we use Holder's inequality and the boundary Poincare 
inequality (I2,14p to obtain with 1 < p^ < oo and 1/qN = 1 — ^-/(^Pn), that 



1 / f%v a da 
P 2 Jn 



< 



< 



1 

P 2 \Jdnndn 2p (x) 



C 



1/PiV 



\Sn? n da 



anndn 2p (x) 
i/pjv 



\u - u x ,2 P \ Pn dy 



p l/p N +2/q N yj mndQ2p{x) 

and now (13.91) follows. 



1/iVp do- 
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Remark. The argument above also gives us a Caccioppoli inequality for solu- 
tions. Let u be a solution of (I2.8P with / = and Jn € L 2 (N), fix x and p and 
let v be as in (|3.7p . Since f vanishes outside f22 P (x), we may use (|2.14p and 
Holder's inequality to obtain J dn \v\ 2 da < Cp f n \\7v\ 2 dy. If dist(x,<9r2) > 2p, 
the boundary term in (I3.7P is zero. When dist(x,(9f2) < 2p, we estimate the 
boundary term in (|3.7p using the Cauchy inequality with e's and the above 
observation to obtain 

I / f N v a da\<C(^ [ \f N \ 2 da + E) + e [ r] 4 \Vu\ 2 dy. 
Jn e JA 2 p(x) Jn 

Using this estimate in (I3.7p and the Cauchy inequality with e's gives 

/ \Vu\ 2 dy<C{\f \u - u x , 2p \ 2 dy + p [ f 2 N da). (3.10) 

JQp(x) P Jn 2p (x) JdQndn 2p (x) 



Corollary 3.11 Let O, D, and L be as in Theorem \3.4\ Let u be a solution of 
the weak mixed problem \2. 8\) with f = on Q,2p(%) and Jn € L 2 (d£l). There 
is a constant C and 70 > so that for 7 in (0,70), we have 

\u{y)-u{z)\<&y^dj \u- 

(3.12) 





f N | 2 da 

(\ 1/2 / \ 1/2 

/ \u\ 2 dy\ +[p[ \fN\ 2 da) ) (3.13) 

Jn 2 p{x) J \ Jdundn 2p (x) J 

The constant 70 depends only on M and C depends on M and 7. 

Proof. We let it be a solution of (j2.8|) with / = in Q,2 P {x). We observe that 
there is a constant c so that if y G ^3 P /2 (x), < s < cp, then ^2s(y) C ^2p(a0 
and we have I a ,2sf = 0. Thus if O s (y) C fJ p (x) and s < cp, the main estimate 
of Theorem 13.41 and (I3.ip imply that for i < min(4,to) 

/ \ 1/* / \ 1/2 / \ 1/2 

4 |V^dy <C((/ |V U | 2 dy +£(/ I/jvI 2 ^ ). 

yn s (y) / vn 2s (2/) / s \JaQnan 2p (x) / 

(3.14) 

Using (|3.14p with the Caccioppoli inequality (|3.10p and Morrey's inequality 
(|2.13p . we obtain the Holder estimate (|3.12p with 70 = min(l — 2/to, 1/2) where 
to is as in Theorem 13.41 

To obtain the estimate (|3.13p . we write 

\u(x)\ < \u(x) - u XiP \ + \u x>p \, 

use the Holder estimate (|3.12p to bound the first term on the right and then 
the Holder inequality. Note that this works even when u x , p = since in this 
case, u vanishes on a nearby piece of the boundary. I 
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We observe a version of the Green identity for solutions of the weak mixed 
problem (I2.8P and the corresponding problem for the adjoint operator, L*. If 
u satisfies Lu = f with du/dv = /n and w is a solution of the adjoint problem 
L*w = g with dw/dv = gw and suppose that / and g are in L P (Q) and /at and 
gN are in L p (dil) for some p > 1. Then, we have 

f u a g a dy- [ u a g%da= [ w a f a dy- [ w a f%da. (3.15) 
Jn Jn Jn Jn 

The identity (I3,15j) follows easily since both sides of (|3.15p are equal to A(u, v). 

The following Corollary is only valid when For the Neumann prob- 

lem, there are non-constant solutions to the homogeneous problem and hence 
the estimates (|3,17p and (|3. 18[) cannot hold. We will give a version of this result 
for the Neumann problem in section [H 

Corollary 3.16 Suppose that D is a non-empty subset of dQ that satisfies 
the corkscrew condition. Let L, $7, and to be as in Theorem \3.4\ suppose that 
t € (2, to), define r by 1/r = 1/2 + 1/t and let r^ = t/2. If u is a weak solution 
of the mixed problem \2. 8\) with /at = and f G L r (£l), then we have 

l|V«|Ut (n) < C\\f\\ L r (n) . (3.17) 

If u is a solution of \2. 8\) with f G W^ 1,1 (fi) and /jv = 0, then we have the 
estimate 

IMI^(n) < (3.18) 

The constant C depends on M, t, and fL 

Proof. Using the estimate (|2.2p . it follows that the map v — > X(v) = Jq f a v a dy 
is an element in W^ 1,2 (Q) with the bound ||A||^-i,2^ < Cr 2 J r ||/||L' (n)- As a 
consequence, if u is a weak solution of the mixed problem (|2.8p with / € L r ({l), 
r > 1 and /jv = 0, we have the energy estimate 

||Vn|| z2(n) <CV2 /r '||/|| ir(n) . (3.19) 

We begin with the main estimate of Theorem 13.41 on domains fl ro (x). We 
choose p with p € (1, r) and apply the estimate (13. ID to obtain 

If iVufdy] <C{rf- l [! \Vu\ 2 dy\ + \\f\\ L r {n) ). 
\Jn ro (x) J \Jn 2ro (x) J 

We use the energy estimate (I3.19P and observe that 1 — 2/t = 2/r'. Finally, 
we cover f2 with domains Q ro (x) and use Minkowski's inequality to obtain the 
first estimate (|3. 17[) . 

The second estimate ()3.18p for solutions of L follows by duality from (I3.15P 
and the first estimate (|3.17p for solutions of the adjoint operator L*. If A* is 
the form for L*, then A*(u,v) = A(v,u), thus it is clear that L* satisfies the 
hypotheses of this theorem exactly when L does. I 
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Our next estimates are valid only when D ^ 0. The estimates may fail for 
the Neumann problem due to non-uniqueness. 

Corollary 3.20 Suppose that D ^ and satisfies the corkscrew condition. If 
u is a solution of the weak mixed problem i2.8\) with f an atom and /jy = 0, 
then we have 

IMIl°°(«) < c. 

Proof. We let u be a solution of the weak mixed problem with / = a, an atom 
for the Hardy space Hp(Q,), and f^ = 0. We suppose that a is supported in 

To estimate u, we fix z in f2, t > 2, and use the Morrey inequality (I2.13P 
and the Holder inequality to obtain 



\u(z)\ < \u(z) - u Z)P \ + \u ZjP \ 



u 



+ P 2/t 1 \\ u \\Lr'(ni< 



(3.21) 



As in Corollary 13. 161 r is defined by \/r = 1/2 + 1/i and r' is the dual exponent 
given by l/r' = 1/2 - 1/t. 

We restrict t to lie in the interval (2, to) with to as in Theorem 13.41 and show 
how to use Corollary 13.161 to estimate the two terms on the right-hand side of 
(|3.2ip . First, we use (|3.17p and the normalization of the atom to conclude that 



||V«|| Lt( n) < C\\a\\ Lr{n) < Cp 1 ' 1 -'. (3.22) 
Next, we claim that for t > 2, we have 

IHItfr-^n) < CP 1 " 27 *- (3.23) 
Given the claim, the estimate (]3. 18|) in Corollary 13.161 implies that we have 

^- 2 l l . (3.24) 



u 



\L T '(Sl) - C P 



To establish the claim (13.23[) . we use that either v X}P = or f^ady = 0, the 
normalization of the atom, and the Morrey inequality f)2. 13j) to write 



f a a v a dy 




/ a a (v - v Xt p) a dy 


< Cp 1 


Jn 




Jn 





which gives the claim (|3.23p . 

The estimate of the Corollary follows from (pmj) . (^22|) . and ([334]) . 



4 The Green function for the mixed prob- 
lem 

In this section, we prove the existence of the Green function for the mixed 
problem and give its main properties. We allow the case when D = dVt which 
gives the Dirichlet problem. 
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We begin by giving our definition of a Green function. This formulation is 
modeled on the definition in Littman, Stampacchia, and Weinberger [18J. We 
say that G(x,y) = (G al3 (x, y)) a ,p=i,...,m ls a Green function with pole at x for 
the mixed problem \2. 8\) if G(x, •) € and whenever u is the weak solution 

of the mixed problem with / G C(O) and /jv = 0, then we have 

u «(x)= [ G^(x,y)f(y)dy. 
Jn 

If we have existence and uniqueness for the mixed problem (|2.8p when the 
right-hand side, / is in C(£l), then it is obvious that the Green function is 
unique. The following theorem shows that the Green function exists and gives 
additional regularity of the Green function. 

Theorem 4.1 Let Q be a bounded Lipschitz domain in the plane and let D C 
d£l be a non-empty open set satisfying the interior corkscrew condition. If L is 
an elliptic operator satisfying \2.1$ , \2. 6\) . and \2. 7\ ) on Wp 2 (ft). Then there 
exists a unique Green function G(x, y) for the mixed problem which satisfies 

G(x,-) £ BMO D (tt), i£!l (4.2) 
V y G(x,-) £ L 2 (n\n r (x)), forxen,re(0,r ) (4.3) 
\G{x,y)\<C{l + \og(d/\x-y\)), x,yefl (4.4) 

\G(x,y)-G(x,z)\<Cf- — %, x,y,z£n, \y - z\ < \\x - #.5) 

\x — y\i 2 

\G(x,y)\ <C dl * t{y > D ^\ x,yen, dist(y, D) < \\x - y\. (4.6) 
\x — y\< 2 

In these estimates d is the diameter o/O and the constants depend on the global 
character of Q as well as M . 

If we let G be the Green function for the mixed problem for L and G the 
Green function for L* then we have 

G^(x,y) = & a (y,x), a,/3 = 1, . . . ,m. (4.7) 

Furthermore, if f € L p (£l) and /jv G L P {N) for somep > 1, then the unique 
solution of the weak mixed problem \2. 8\) is given by 

u«(x)= [ G Q P(x,y)f(y)dy- [ G a ^x,y)f N (y) da. (4.8) 

Proof. We fix x € ft, a in {1, ...,m}, and p E (0,ro). We let G^'(x, ■) be 
the solution of mixed problem (j2.9j) for L* with this g = e a Xn p (x)/\Qp( x )\ an d 
9N = 0. We fix / in L p {9.) and f N G L PN (N) with p > 1 and^iv > 1 and let 
u be a weak solution of the mixed problem for L (|2.8p with data / and f^- If 
we let w = Gp'(x, •) in (13.15j) . we obtain 

/ u a (y)dy= [ Gf(x,y)f(y)dy- [ Gf (x,y)f N (y) da. (4.9) 
Jn p (x) Jn Jn 
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If /jv is zero and / is an atom for Hjj(Q), we may use the estimate of Corollary 
13.201 and the equivalent norm on BM0e>(£1) given in (|2.10p to conclude that 
-)||*,d < C with the constant C depending only on M and the global 
properties of Vt. The Banach-Alaoglu theorem gives that for each x, there is 
a sequence {pj} with lim^oo pj = and a function G(x, •) G BMOd(Q) so 
that Gp.(x, •) converges to G(x, •) in the weak-* topology of BMOd(SI). Since 
u is Holder continuous, the left-hand side of (|4.9p converges to u a (x). Since 
L P (Q) C Hjj(fl) for all p > 1, we obtain the representation (|4.8p in the case 
that /n = 0. This gives us that G(x, •) is a Green function for the mixed 
problem with pole at x. Thus, we have established that the Green function lies 
in BMO D (n), (142]) . 

If we choose any sequence {pu} with lim^oo p^ = 0, the above argument, 
applied to the rows of G Pk , gives a subsequence of {G Pk (x, •)} which converges 
to a Green function. As the Green function is unique, the limit must be the 
function G(x,-). This implies that the entire family {G p (x, -)} p converges to 
G(x, •) in the weak-* topology of BMOd{^)- 

Next, we recall that BMOd{Q) C L p (VL) for any p < oo. Thus we may use 
the Caccioppoli inequality (I3.10|) . to conclude that 

/ \V y G p (x,y)\ 2 dy < C(r), p < r/2. 
Jn\n r (x) 

This estimate will also hold for the limit and thus we obtain the conclusion fj4.3|) 
and that the rows of G(x, ■) are solutions of L*G a '(x, •) = in £1 \ {x}. More 
precisely, we have A(<j), G a ' (x, •)) = whenever cj) 6 W^ 2 {VL) and vanishes 
in a neighborhood of x. Since G(x, •) is a solution of the adjoint equation, we 
have the estimates of Corollary 13.111 in \ {x}. 

We show how to use these estimates to obtain the pointwise bounds of the 
Theorem. If x and y are in f2, we may find a chain of domains flj = Vt Pj (yj) for 
j = 0, . . . ,N so that a) % C %n% +i 7^ for j = 0, . . . , N-l, b) p > 

c\x - y\, y = y, p N = r /2, c) 1 < Pj+i/pj < 2, and d) N < Clog(d/\x - y\). 
Since G(x, •) is in BMOd(£1), we have that 



\*,D- 



I G(x,y)dy--f G(x,y)dy <C\\G(x,- 

Since the rows of G(x, •) are solutions of the equation L*u = 0, the properties 
of the chain {%} and the bound (I3.13P implies the pointwise bound (I4.4p for 
G. To obtain the Holder continuity (|4.5p . we use that G lies in BMOd({1) and 
the local Holder estimate (13. 12f) . If we fix x and y, then we may apply the local 
Holder estimate on a local domain £l p {y) with radius p comparable to \x — y\. 
The boundary estimate (14. 6p follows immediately from (14. 5p . 

Next, we turn to the symmetry property of the Green function (|4.7p . We 
let G p and G p be the approximate Green functions for L and L* used in the 
construction of the Green function. Using the Green identity (|3. 15j) we obtain 

/ Gf(x,z)dz = -f Gf(y,z)dz. 
Jn p ( y ) Jn p (x) 
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Using the Holder continuity of G p and G p in the second variable and the Arzela- 
Ascoli theorem, we may extract a subsequence which converges uniformly on 
compact subsets of Q \ {x}. Letting p tend to zero, we obtain G a ^(x,y) = 
& a (y,x). 

Finally, to obtain the representation formula of the Theorem for solutions 
with / and Jm not zero, we may use the Arzela-Ascoli theorem to find a sequence 
G p . (x, •) which converges uniformly on dQ. Thus, we may take the limit in (]4.9j) 
to obtain the representation formula of the Theorem. The convergence of the 
integral on f2 follows since L P (Q) C Hjj(Q.) and G p (x,-) converges weakly in 
BMO D {n). ' I 



5 The Green function for the Neumann prob- 
lem 

In this section we consider the Green function for the Neumann problem (which 
is the mixed problem in the extreme case where D = 0). Most of our arguments 
parallel the construction of the Green function for the mixed problem. However, 
there is an additional complication. The homogeneous Neumann problem for L 
has non-trivial solutions. Hence, we need to impose compatibility conditions on 
the data and conditions to guarantee uniqueness of solutions. It seems that the 
most natural condition for uniqueness involves the boundary values of solutions. 

We let L be an operator as defined in (|2.3p and we consider the form A 
now defined on W^' 2 (Q), the homogeneous Sobolev space of functions with 
one derivative in L 2 (Q). Since we have chosen to norm this space by the 
expression (J*^ \Vu\ 2 dy) 1 / 2 , the elements of this space will be equivalence classes 
of functions under the equivalence relation u is equivalent to v if u — v is a 
constant. 

We let V = {v: A{v,<j>) = for all <j> G W^ 2 (fl)} and V* = {v : A(<j>,v) = 
for all 4> G Wq' (tt)} denote the solut ions of the homogeneous Neumann prob- 
lems for L and L* respectively. Under our ellipticity assumption (|2.6p and 
boundedness of the coefficients (12.41) we have that these spaces are finite di- 
mensional and the Fredholm alternative implies that dimV = dimV*. The 
Holder estimate of Corollary 13.111 implies that the elements of these spaces are 
Holder continuous. 

We give the weak formulation of the Neumann problem for L 

A(u, 4>) = -(f, 0) + (f N , 4>) dQ , 4> g w^ 2 (n) 

u G W^ 2 (n) (5.1) 
fan uOLva da = °> v eV. 

If v lies in V*, then we have A(u,v) = for u G W^' 2 (Q). Thus, if we are to 
find a solution to (|5.ip . we must have that / and /n satisfy the compatibility 
condition 

(fN,v)aa-(f,v)=0, veV*. (5.2) 
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For the operators we consider, the constant functions lie in V*. Thus, the first 
line of (15. ip is satisfied if is any representative of an equivalence class in 

We will use the following technical result when we construct the Green 
function. 

Proposition 5.3 Suppose Q is a Lipschitz domain and that the operator L 
satisfies \2.1$ , \2. 6\) . and the coerciveness condition (2.1) with D = d£l. 



a) The norms |[ • ||^2(q) + || • H^ 1 , 2 ^) and \\ ■ \\L 2 (dQ.) are equivalent on the 
space V . 

b) If (i = (fj, , . . . , fi m ) is a finite R m -valued Borel measure on Cl, then we 
may find A^ 6 V so that 



Jn Jan 



v a d\i a = / \y a da, v £ V. 



We have the estimate \\Xn\\L 2 (dO.) — ^IImII where \\n\\ denotes the total variation 
of ii. 

Proof. Since V C W ' (fl), it follows that f m \v\ 2 da is finite on V. If v e V 
and J S q \ v 1 2 da = 0, then v = on dtl and it follows that v is in the Sobolev 
space Wg^(Q). Since we assume that the form A is coercive on this space, it 
follows that v = 0. Thus, we have || ■ ||z / 2 (an) 1S a norm on this space. Since 
V is finite dimensional, it follows that the norms || • ||l 2 (q) + || • H^ 1 * 2 ^ and 

II ' lli 2 (<9C) are equivalent. 

To establish part b), observe that the local boundedness estimate (I3.13P 
implies that v — > Jjjti d[i a is a continuous linear functional on V. According 
to part a), V is a Hilbert space under the inner product f dn u a v a da. Hence, 
the Riesz representation implies that we have a unique A^ € V so that 

v a d f i a = [ \ a v a da. 
n J do, 

Corollary 13.111 implies that the elements of V are bounded functions, the esti- 
mate for A^ follows from the Riesz representation theorem. I 

Next, we give a standard existence theorem for the weak Neumann problem 

4SU). 

Theorem 5.4 Let fi be a Lipschitz domain and suppose that L satisfies the 
ellipticity condition \2. 6\) . has bounded coefficients {2.1$ , and that {2. 7\ ) holds 
for D = dO,. If f and f^ satisfy the compatibility condition h5.2}) . then the 
weak Neumann problem h5.1\) has a unique solution and the solution u satisfies 
the estimate 



Wh^n) + l|u|lifc}.a (n) - c (Wf\\w- 1 - 2 {n) + WfN\\ w -i/2,2 



(dSl)' 
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Proof. The existence of solutions is a standard consequence of the Fredholm 
alternative. See Gilbarg and Trudinger |12} Theorem 8.6] for the proof of 
a similar result. To establish uniqueness, if u solves (|5.ip with / = and 
In = 0, then we have that it € V. Part a) of Lemma 15.31 and the condition 
Idfi uOLva do~ = for v G V imply that u = 0. I 

We now give estimates for solutions of the Neumann problem to take the 
place of Corollary 13.161 Since the homogeneous Neumann problem may have 
non-constant solutions, Corollary 13. 161 cannot hold for all solutions of the Neu- 
mann problem. 

We consider a function / in L r (Q), r > 1, and use Lemma 15.31 to find 
a function A/ G V* so that the pair / and \f\do, satisfy the compatibility 
condition needed to solve the Neumann problem, 

(X f ,v) du -(f,v)=0, vGV*. (5.5) 

We let A f be defined by A/(0) = f m \]4> a da - J Q f a (j) a dy. We observe that 
the constant functions lie in V* and hence we have Af(c) = for all constant 
vectors c G R m . The elements in V* are bounded and for all p, 1 < p < oo, we 
have the inequalities 

r o 2/P W u ~ u\\lp(Q) + - u\\ LP{dn) < C\\Vu\\ L 2 {n)j (5.6) 

where u = f^udy. It follows that Af lies in W~ 1,2 (Q) and thus the solution of 
the Neumann problem 

A(u,cp)=A f (cp), c^eW^ip) 
u G W^ 2 (n) (5.7) 
Ion u0lya do- = 0, v £ V. 

exists and satisfies ||Vu||x2m) < C(r,ro)\\f\\L r (a,)- Next, we consider the adjoint 
problem for g G L r (f2), 

' A(<t>,w)=A*(cj>), <p ^ W l ' 2 {Q) 
< weW^' 2 (Q) (5.8) 
k J m w a v a da = 0, v G V*. 

where A*((p) = (X g , <p)sn - (g, 4>) and \ G V satisfies J Q v a g a dy = f m v a Xg da, 
v G V. Similar considerations give the existence of a solution to the adjoint 
problem with ||Vw;||£2(fi) < C||s r ||i,*-(f2) - Furthermore, from (I3.15|) we obtain 



u a g a dy= / w a f a dy. (5.9) 
n Jn 

The boundary integral Jq^XJuj" da vanishes since the Neumann data Xf lies 
in V* and w is perpendicular to this space in L 2 (dQ). The other boundary 
integral vanishes for similar reasons. 
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Corollary 5.10 Let t G (2, to) with to as in Theorem \3.4\ and r be defined by 
1/r = 1/2 + 1/t. If f G L r (fl), then the solution of the weak Neumann problem 
|5. 7| ) satisfies the estimates 

||V«|| L t (0 ) < C||/|Ur (n) (5.11) 

and 

h\\ L r> {n) < C\\f\\ w ^ t , {n) . (5.12) 

Proof. With the work above the proof is the same as Corollary 13.161 The 
proof of Corollary 13.161 fails at the first line, because we do not have a Poincare 
inequality. Since Af (u) = A f {u — u) , we may use the estimates (|5.6I) to show 
that A f lies in the dual of W^' 2 (Cl). With Af G W^ ,2 (Q), the proof of this 
Corollary is identical to the proof of Corollary 13.161 I 



Remark. For the second estimate of Corollary 15. 101 to be useful, we must have 
that / lies in W~ 1,1 (CI). In particular, we must have that the mean- value of / 
is zero. 

Corollary 5.13 If a is an atom for the Hardy space H^(fl), then the solution 
of J 5. 7| ) with f = a satisfies 



Mlz,°°(n) 



< C. 



Proof. Given the estimates Corollary I5.10| the result follows as in Corollary 
I3T201 I 

We are ready to define a Green function for the the Neumann problem. We 
say that G(x, ■) is a Green function for the Neumann problem with pole at x if 
G(x, ■) is in L 1 (fi) and whenever a is an atom and u the corresponding solution 
to (|5.7p with / = a, then we have 

u a (x)= [ G a P(x,y)aP(y)dy. 

It is clear that the Green function is unique up to a constant. 

Theorem 5.14 Let CI be a Lipschitz domain and suppose that L satisfies the 
ellipticity condition 112.6}) . has bounded coefficients {2.1$, and the form for L is 



12 

coercive on Wqq(Q). Then there exists a unique Green function G(x,y) which 
satisfies the following estimates 

G(x,-) G BMO$(Cl) (5.15) 
V y G(x,-) £ L 2 (Cl\n r (x)), forxeCl,r>0 (5.16) 
\G(x,y)\<C(l + \og(d/\x-y\)), x,y£Cl (5.17) 

\y — £ — l 

\G(x,y) - G(x,z)\ <C _ , x,y,zefl, \y - z\ < ~\x - y[5.18) 

\x y\ i 
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If G and G are the Green functions for L and L* , respectively, we may find 
representatives of G and G so that 

G a ?{x,y) = & a (y,x). (5.19) 

Furthermore, if u is a solution of the weak Neumann problem with f £ LP(£l) 
and /at € L p (dQ) for some p > 1 and satisfying the compatibility condtion 
if 5.2i). then we have 



u °(x)= / G^{x,y)f{y)dy- / G a ?(x, y)f N (y) da, a = l,...,m. 
Jn Jan 

(5.20) 

Proof. We fix x £ n, a £ {1, ...,m}, p with < p < rp, and let g = 
e aXn p (x)/\Qp( x )\ where e a is the unit vector in the direction of the crth coor- 
dinate axis. We let A"' € V be chosen so that g and A"'|an satisfy the com- 
patibility condition needed to solve the adjoint problem (|5.8p . We let G p '(x, ■) 
be the solution of (15.81) with g as above and gpj replaced by A"'. We let u be 
a solution of the Neumann problem f)5. If) with data / and /jv in L p (£l) and 
L p (dQ), p > 1, respectively. From (I3.15p . we have 

/ u a (y)dy= [ Gf(x,y)f^(y)dy- [ Gf(x,y)a^y)dy (5.21) 
Ju p {x) Jdn Jn 

If we let / be an atom and /at = 0, then Corollary 15.131 implies the left-hand 
side of (I5.2ip is bounded by a constant that is independent of a. It follows that 
||Crp"(a;, ')ll*,0 < C with C independent of p. We may use compactness in the 
finite dimensional space V and the Banach-Alaoglu theorem to find A Q ' £ V, 
G a '{x,-) € SMO0(O), and a sequence {/Oj}"^ with linij^oo p.,- = so that 
Gp' (x,-) converges weakly to G a '(x, •) in BMO${Q) and A£\ converges to A"' 
in L 2 (dQ). Combining this weak convergence of G Pj in BMOq(CI) with the 
Holder continuity of u, we obtain that 

u a {x)= [ G&{x,y)aP{y)dy. 
Jn 

Thus, we have found our Green function and we have (|5.15p . From the unique- 
ness for the solution u, it follows that the limit G(x, ■) is unique and thus we 
have that G p converges for all p and not just a subsequence. 
Using the Caccioppoli inequality ()3.10p we can show that 

/ \VG p (x,y)\ 2 dy <C(r), p<r/2. 
Jn\n r (x) 

As this estimate is uniform in p for p sufficiently small, we obtain that for all 
4> € Wq ,2 {VL) which vanish in neighborhood of x that 
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The pointwise estimates (|5.17p and (|5.18p follow from (|5. 15[) and the local 
estimates in Corollary 13.111 As G has non-zero Neumann data, it is important 
that these results allow non-zero Neumann data. The argument is identical to 
that of Theorem 14.11 

We establish the symmetry property (I5.19|) . We let G p and G p be approx- 
imate Green functions for L and L* as defined above. From (|5.9p . we obtain 
that 

/ G p la (y,z)dz = -f Gf(x,z)dz. 
Jn p (x) J^p(y) 

Now, we may let p tend to zero and obtain (|5.19p . 

Next, we claim that the mean value f^G p ^(x,y)dy is bounded for all p. 
As a first step, let u be the solution of (15.71) with / = epxn- According to 
Corollary I5.10| the solution u lies in L r (£l) for some r. As the Neumann data 
Xf lies in V* and hence is bounded, we may use the estimates of Theorem 13.41 
to conclude that Vu lies in L for some t > 2. Since u lies in Z/(f2) and Vu 
lies in L (O), we may conclude that u is bounded. We apply (|3.15p and obtain 

/ u a dy= [ Gf(x,y)dy. 
Ja. p {x) Jn 

Since u is bounded, the claim follows. Since G p '(x, •) is bounded in BMO^(Q) 
and the mean values are bounded, it follows that a subsequence of G p '(x-) 
converges weakly in L P (Q) for each p finite. In addition, G p '(x, •) is a solution 
of L*w = in a neighborhood of the boundary and the Neumann data \ p lies 
in V and hence is bounded. Thus, we may extract a subsequence {G Pj (x,-)} 
which converges uniformly on d£l. We may let p — > + in the representation 
formula (|5.2ip and use the continuity of u to obtain the representation formula 



6 The Green function in the plane 

In this section, we define a Green function in the plane. We will work in 
the homogeneous Sobolev space W /1 ' 2 (R 2 ) which consists of functions <p with 
V(/> € L 2 (R 2 ). We norm this space with 

IMIwn, 2( R 2) = (J Wfdy} ' • 

and the elements of this Hilbert space will be equivalence classes of functions 
under the relation u is equivalent to v if u — v is constant. In R 2 , the local 
domains Q p (x) are disks and we will use the more standard B p {x) = {y : 
\x — y\ < p} to denote these disks. 

We assume that the form A is coercive in the sense that 

A(u, u) > M _1 / \Vu\ 2 dy, u G W 1,2 (R 2 ). (6.1) 
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We may see that this condition holds for the Lame system using an approxi- 
mation argument and integration by parts. 

It is an immediate consequence of the Lax-Milgram theorem that the weak 
formulation of the problem in the plane 

■ A(u, # = ^e^ 2 (R 2 ) > 

u £ W 1 ' 2 (B?) { ' 

has a unique solution when / is in the dual of H^ 1,2 (R 2 ), W~ 1,2 (R?). 

We will approximate u by considering the Dirichlet problem in disks Br = 
{x : \x\ < R} for R > 0. We let a be an atom for R 2 that is supported in B p {x) 
and for R large, we let ur be the solution of the Dirichlet problem 

f A(u R , <t>) = -J n dy 4> € (B R ) 

{ ur£W^ r (Br) 

From Proposition l3.3l and Corollary [3T201 we have that ||uij||£,oo^ i s ii )+|| ^ u r\\l 2 (b r ) < 
C and from Corollary 13.161 and Morrey's inequality (|2.13p 

W u R\\Lr'(B R ) + SU P | T _ 7 ,i 7 < °(P)- ( 6 - 3 ) 

The Holder index 7 = 1 — 2/t with t and r' as in Corollarv l3.161 The estimates of 
Corollary 13. 16l are scale invariant and thus hold uniformly in R. The dependence 
on p arises because the norm of an atom in U'(Br) and W~ 1,1 (Br) will depend 
on p. Thus, we have a function u £ M^ 1,2 (R 2 ) so that lim/j_ s>00 ur = u weakly 
in W 1,2 (Bs) for each S > and u solves (|6.2p with / = a. Note that since the 
limiting function u is unique, we have convergence for the entire family not just 
a subsequence. Furthermore from the Rellich compactness theorem, we have 
that ur converges in If (B$) for each S > 0. The Holder estimate in (I6.3P and 
the Arzela-Ascoli theorem imply that we also have that ur converges locally 
uniformly to u. According to Corollary 13.201 the functions ur are uniformly 
bounded, hence the same holds for u. Thus, if / = a, an atom, the solution of 
(|6.2p may be chosen so that 

< C (6.4) 
NI^CR^+Bup '^r.y 1 < C(p). (6.5) 



We give a definition of the Green function in the plane. We say that G(x, ■ 
is a Green function in the plane for \6.2i) with pole at x if G(x, ■) is in L ; 1 oc (R 2 ) 
and for each atom a, the solution of (|6.2p is given by 



u a (x)= G al5 (x,y)a l3 (y)dy, a = l,...,m. (6.6) 

Since solutions of this weak problem are unique, it is immediate that for each 
x, G(x, •) is a unique element of BMO(H 2 ). In other words, G(x, ■) is unique 
up to a constant. 
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Finally, we give a theorem which establishes existence and regularity of the 
Green function in R 2 . 

Theorem 6.7 If L satisfies {2.1$ and \6. then there is a unique Green func- 
tion G(x, ■) with pole at x and the Green function satisfies 

G(x, •) G BMO(R 2 ) (6.8) 
V y G(x,-)€L 2 oc (R 2 \{x}) (6.9) 

\G(x,y)--f G(x,z)dz\ <C(l + |log(|ar-y|)|), x,y eB 2 (6.10) 

\G(x,y) - G(x,z)\ < C^^, if \y-z\< \\x - y\ (6.11) 

\x — y\' 2 

The exponent 7 is as in Corollary \3.11\ and C may be chosen to depend on M 
and 7. 

Finally, if G and G are the Green functions in the plane for the operators L 
and L* , then we may find representatives which satisfy the symmetry condition 

G a P(x,y) = & a (y,x). (6.12) 

Proof. To construct the Green function G(x, •) we fix x, p > 0, and let f p = 
^zXb p (x)-§^Xb 2/p (x)\b 1/p (x)- As f p e a is an atom (though the constant depends 
on p), we may let G p (x, ■) be the weak solution of the equation L*G p '(x, •) = 
f p e a . We let a be an atom and u the solution of Lu = a constructed above. 
As in the proof of (|3.15p , the weak formulations of the equations satisfied by u 
and G p '(x, •) give 

A(u,G a p (x r ))= [ Gf(x,y)aP(y)dy=[ f p u a dy. (6.13) 
Jr 2 Jr 2 

The estimate for |M| i0 o( R 2) in (|6.4p implies that | J G p ^ (x, ^a 13 dy\ = \ J f p u a dy\ < 
C . Now, we may conclude from (f2"TT0D that G p (x,-) lies in BMO(R 2 ). Thus, 
we may use the Banach-Alaoglu Theorem to find a function G(x, •) and a se- 
quence {G Pk (x, •)} with lim^oo = so that G Pk (x, ■) converges to G(x, ■) in 
the weak-* topology of BMO(R?). From the estimates in (|6.5p . it follows that 

u a (x) = lim / f p u a dy. 

Hence, we obtain the representation formula (16.6p . We may use Caccioppoli's 
inequality (|3.10p to obtain uniform bounds on S/ y Gn(x,-) in L 2 of compact 
subsets of R 2 \ {x}. Thus G(x, ■) is a solution of L*G(x, •) = in R 2 \ {x}. 
Now the pointwise estimates (|6.1U|) and (|6.11|) follow as they do for the mixed 
problem. 

Finally, we establish the symmetry property (|6.12p . As our construction of 
G, stands we have no information about the behavior of G in the first variable. 
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We begin by claiming that we can fix a representative of G(x, •) so that G is 
locally integrable in R 2 x R 2 . Towards this end, we observe that if we fix x 
in R 2 and let h PjX = ^z(xb p {x) ~ Xb p (o)), then e a h PyX lies in ^"^'(R 2 ), the 
dual of W^CR 2 ), t > 2, with ||/i p e a ||^_i it / (R2) < C^ 1 " 2 / 1 . We let v p (x, •) be 
the solution of L*v p {x,-) = e a h PyX . Using the estimates of Corollary 13.161 we 
may show that the map x — > v p (x, •) is a continuous map from R 2 into L r (R 2 ) 
where r' is as in (13.180 . If we fix a representative of G p (0, •), the function 
G p (0, •) + v p (x, •) gives an approximate Green function that is locally integrable 
in R 2 x R 2 . If we let p — > + , we obtain the same conclusion for G. 

We let G and G be the Green functions as constructed in the previous 
paragraph for L and L*. We fix atoms a and b and let u and v solve the 
equations Lu = a and L*v = b. From the weak formulation (|6.2p . we have 
A(u, v ) = f u a b a dy = f arv^ 1 dy. Using the representation (16. 6j) and Fubini's 
theorem, we obtain 

/ b a (x)G a ^(x,y)a l3 (y)dxdy = [ b a (x)G^ a (y,x)a^(y)dxdy. 

JR 2 xR 2 JR 2 xR 2 

As this holds for all atoms a and b, we have functions 4> an d (j> so that 

G(x,y) + <t>{x) = G(y,x) + cf>(y). 
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